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Abstract 

We introduce a class of nonlinear damping models with application to flexible flight 
structures characterized by low damping. We are able to obtain approximate solutions of 
engineering interest for our model using the classical “averaging” technique of Krylov and 
Bogoliubov. The results should be considered preliminary pending further investigation. 
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1. Introduction 


The problem of characterizing the damping mechanism in flexible structures has 
received renewed attention in recent years in connection with the need to stabilize flexible 
flight structures such as antennas deployed in space. The damping models even when 
simplified to be linear appear to lead to rather complex mathematics if the structure is 
described by partial differential equations and much progress has been made (the analyticity 
of the generated semigroup has been shown to be essential). But experimental evidence as 
in SCOLE [ 1 ] seems to support the need for nonlinear models — the decrement is much 
smaller than predicted by linear models. Some of the difficulty inherent in handling 
nonlinear models is offset by the fact that damping, whatever its nature, is still small. This 
opens up in particular the feasibility of obtaining approximate solutions using the classical 
averaging method of Krylov-Bogoliubov [2], 

In this paper we study a class of nonlinear models and approximate the response by 
the Krylov-Bogoliubov technique. We use a modal expansion and neglect off-diagonal 
terms. The emphasis is on useful engineering solutions rather than abstract mathematics. 

We begin in Section 2 with the primary nonlinear damping model for the simplest 
system — the one-dimenisonal or single-mode case. We emphasize in particular one 
feature that emerges, viz., the potential lack of identifiability from response data. In Section 
3 we generalize to the multi- (non-finite-) dimensional case. In Section 4 we show the 
relevance of the Krylov-Bogoliubov technique for approximating solutions to nonlinear 
boundary feedback. We may mention that there is much work — even classical in nature 
— on nonlinear oscillations such as the nonlinear pendulum where the spring constant is no 
longer linear; however, relatively little attention appears to have been paid to the small 
nonlinear damping term case. 




2. Single-mode Example 


To illustrate ideas, let us begin with a one-dimensional (single-mode) example: 

x(t) + zD(x, x) + CO 2 x(t) = 0 (2.1) 

where the dots indicate time-derivatives, as usual. We assume that: 

D(x, x)x £ 0 (2.2) 

so that for E(t), the energy 

E(t) = j (x(t) 2 + (£> 2 x(t) 2 ) (2.3) 

we have 

J-£(0 = -eD(x,x)x < 0 (2.4) 

satisfying the energy nonincrease requirement. The particular choice for D(x, i) we 
shall make is: 

D(x, x) = 2co£i + yx 2m |x| a iK 2 " + 1 )[xp (2.5) 

where m, n are nonnegative integers, 

0 < a, p and 0 < a + P < 1; 0 < £ < 1 , 0 < y < 1 . 

For small enough e we may apply the averaging method of Krylov-Bogoliubov [2, 5]. 
Thus, we write for the approximate solution: 

x(t) = a(t) sin (cor + 4>(r)) (2.6) 

where the amplitude function a(t) and the phase function <|>(r) are slowly varying over 
the period T = 2jt/co. According to the K-B approximation [2]: 
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Now we can readily calculate that for our choice, because of (2.2), 
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and is a constant less than 1/2. Hence letting p = 2m + 2n, we have 

z?+l + a+B 2n + B ^ 

Ji - -e(co£a + <r P G) p y|i) 

We may set e = 1 without loss of generality since we may absorb it into £ 
Then 
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yielding 


a(t) = a(0)e'“ ; [l + a(0) p+ “ + P co 2 "" 1 + P ^ (1 - e _, “ C(,,+ a + P) )] 


’ p + a 4 P 


We can readily verify that for £ = 0, we have 


a(t) = a(0)[l + a(0) />+a + * 3 co 2n + PY|i(/> + a + (3)r] f + a + P 
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( 2 . 10 ) 


d( j) 
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and y. 

(2.13) 

(2.14) 


( 2 . 15 ) 



The case y = 0 is even more obvious. One salient fact that emerges immediately from 
(2.14) and (2.15) is that it would be difficult to resolve p + a + (3 into its components 
from response data, unless we can change co. 

Note also from (2.15) that the rate of decay is not exponential in t and further the 
decrement over any integral multiple of the period depends on the initial amplitude as well 
as the frequency of oscillation. Finally for inegral k and 

_ 
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we have, taking logarithms and setting c = p + a + (3 
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For small C, this is well approximated by 

-2nkC, - ^ log (1 + 2nkXa(0) c ) , 


where 


X = 


2 n + p - 1 

CO p 


Ylt . 


The slope (as a function of k) 


r 27ciU(0) c 

c(l + 2KkXa(0) c ) 


(2.17) 


and hence the linear damping term is yielded by the asymptotic slope as £-*<*>, while 
for small k there is a marked curvature which depends also on the initial amplitude 
a(0). The inital (at k = 0) slope 

= -2tc£ - 27t ^°) c - (2.18) 

is larger (in absolute value). The second derivative being positive, the curve is convex 
— CUP. This is in excellent qualitative agreement with SCOLE damping data: see 


Figure 1 where amplitude is plotted on logarithmic scale (period = 5 seconds). 


SCOLE DAMPING 
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ORJQJNAL PAGE fS 
OF POOR QUALITY 


To get another version of (2.16) we may replace (2.7) by the more exact formula 
qit+T} - a(t) _ ~^K 0 (a(t)) (2.19) 

and hence using 
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3. Multidimensional Generalization 

Analogous entirely to the one-dimensional case, we may write the general nonlinear 
dynamic equation for flexible structures [2] as 

Mx{t) + D(x(t), x(t)) + Ax{t) = 0 (3.1) 

where the state x{t) ranges in a separable (real) Hilbert space H; M is a self- 
adjoint positive definite (with bounded inverse) operator on X onto X; A is a self- 
adjoint nonnegative definite closed linear operator with domain dense in X and with 
compact resolvent; we shall (for simplicity) assume that zero is in the resolvent set of 
A. In the linear case 

$(*(0> *(0) = Dx(t) (3.2) 

where D is also a self-adjoint nonnegative definite closed linear operator whose 
domain includes that of V/4 . In the most important case we further specify that 
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D$ k = (3.3) 

where {<}>* } are the M-orthonormalized eigenfunctions of A with eigenvalues co* such 
that 

Aty k = <o 2 k M$ k . (3.4) 

Here ^ k is the damping ratio. If = £ and we have strict proportional damping — 
see [3] for more — D is then essentially the positive square root of A (except for M). 
More generally we require that lim £ £ > 0. In the nonlinear analogue of (2.5) we 

k 

set 

3X4»>, ♦*) = 0 j* k (3.5) 


and more generally for x, y such that 



2 ,2 4m + 2 a , 2 n + 2 + 2 P 2 

I to* b k + I a k b k y k < oo 

(3.6) 

where 

a k = l x > <$>*] ; b k = [y, 4>*] 

(3.7) 

we define: 

[%(x,y), = y k al m \af b 2 k n+l \b k f + 2 ^( 0 k b k 

(3.8) 

where, as before, 

m and n are nonnegative integers and that 



0 < a , P<1; a + P < 1 ; 0<y. 

(3.9) 

Note that 

[2)(x, y), y) > 0 


for every x and 

y. Hence if 



E{t) = j{[Ax(t), j:(r)] + [Mx(t), i(r)]} 


we have that 

J-£(r) = -[3)(x(0, *(')), i(/))] < 0 . 

(3.10) 


Or, the energy is nonincreasing. Using the modal expansion 

a 
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x(t) = X a k (t) <(>* (3.11) 

we see that for each k 

d k {t) + 0 i 2 k a k (t) + a k (t) 2 \a k (t)\ a a k (t) |a*(f)| P + 2^ k (0 k d k (t) = 0. (3.12) 

We can therefore invoke the K-B averaging procedure obtaining the approximate solution 

a k (t) = A k (t) sin (co*f + <|>*) (3.13) 

A k (t) = A k (0) e"' ; *“*[l + A*(0) c (D*' , + P-1 ^(1 - e"' 5 *®**)] ‘ . (3.14) 


And for t, k = 0, 

A k (t) = A*(0)[l + A k (0) c «> 2 k n + t Y*jtcr] ‘ 

where, as before, 

. 2* Q 

jj, = — J sin 2m 4> cos 2rt <t> |sin <|>| a looser cos 2 <(> d<() . (3.15) 

o 

For a + p = 0, we can give a kernel representation. Thus 

Z = 3)(x, y) = 1 Y, $,[<!>„ xf m y] 2 " + 1 

where 

I yf < 00 ; y, ^ 0 ; 

and for the concrete realization K = £. 2 ( 0 , L), the corresponding “kernel” would be 
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which will satisfy (3.10), since the y, in (3.16) are nonnegative. In the notation of 
(3), the “roll” equations for example will have the form: 


pA u^it, s) + ElqUp'Xt, s) 




d 3 u <) (t, s) 
dtps 


- J ' J W(s,a 


l > 
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dtd 2 S 


d ° i dG 2m 

= 0 . 


ds\ •" ds 2n + i 


It is clear that we may generalize (3.17) without recourse to modes. The “nonlocal” 
nature of the operator should hardly be surprising, since this is already so in the 
linear case if we want strict proportionality (£* « £) for example. 


4. Application to Nonlinear Boundary Feedback 

In this section we shall apply the K-B averaging technique to obtain approximate 
solution to the response of a flexible structure to nonlinear boundary feedback control. The 
control effort is small so that the K-B approximation is reasonable. We follow [4] for the 
model where the “boundary” is finite-dimensional. Thus we have in the same setting as 
Section 3, but omitting the natural damping term: 

Mx{t) + Bf(B*x(t)) + Ax(t) = 0 (4.1) 

where B means R m onto X and /(•) maps R m into R m and is such that 

[/(«), u] > 0 for u * 0 . (4.2) 

Using the modal expansion as in Section 3: 

x(t) = Ia*(r)<t>* 

we obtain 

d k {t) + CO* a k {t) + [/(X aj(t)bj), b k ] = 0 (4.3) 
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where 


B**k - t>k ■ 


Taking the approximation 


[/(I aj(t)bj), b k ] = [f(a k (t)b k ), b k ] 


we see that setting 


where 


a k (t ) = A k (t) sin (co*r + <^*( 0 ) 


£♦*«> - 0 

i . ... *>frM 0 ) 

5T' 4 * w ’ - o>. 


- 2 n 

AT°(a) = 2jr J cos <t>h t ), h*] cos <|> d§ 


To simplify matters further let us take 


m - 1 . 


We shall take: 


[f(avs k b k cos <t>), b k ] = b k f(ab k { 0 t cos <j>) 


/(«) = A. tan - 1 u 


which is consistent with ( 4 . 2 ). Then ( 4 . 3 ) becomes 


Hence 


ci k ( t) + to t a*(r) + ^ ^(0 — 0 


1 2,1 

/C 0 (a) = J tan -1 {ab k 0 ) t cos <|>) cos (j) <f<|) 


= + aZb Wk ~ !] 


A k (t) dA k (t) 
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To solve this, let 


z(r) = 1 + A k {t?bl<i>l 
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so that 


-7P r = -b\Xdt 

- 1 ) 


(4.6) 


Let 


F(z) = e 1 (V7 - 1) , z > 1 . 


(4.7) 


Then 


F'(z) > 0 


for z > 1 


and hence we may define the inverse function 


F(z) = y ; 2 = F 1 (y) . 


Thus (4.6) has the solution: 


2(r) = F _1 [F(2(0))e' fc * X '] , 


(4.8) 


b k (O k A k (t) = Vf-’(F(1 + A k {0) 2 bW k )e~ blkX ') - 1 (4.9) 


where 


z(0) = 1 + A k (0) 2 b\ G>1 > 1 


unless (0) = 0. Note that 


F (y) -» 1 as y -» 0 


and hence z(t) decreases monotonically to 1 and hence the amplitude A k (t) decays to 
zero asymptotically. 

Note that the decay rate depends on the control effort Xb\ as well as the initial 
amplitude. Of course we have in (4.4) yet another nonlinear damping model. Following 
(2.22) we have: 

lo a j + i _ 2 nX [ VTTggg - 1 ] 

8 a J 2 ( 0 1 


= A k (jT) ; T = 


2_tt 


6S - 


where 
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